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soil moisture from a generalized algorithm requiring only received 
power and the mean elevation of the test site. The results demonstrate 
that the soil moisture of about 90 percent of the 20-m by 20-m pixel 
elements can be predicted with an accuracy of ±20 percent of field 
capacity within relatively flat agricultural portions of the test site. 
Radar resolutions of 93 m by 100 m with 23 looks or coarser gave the 
best results, largely because of the effects of signal fading. For 
the distributions of land-cover categories, soils, and elevation in 
the test site, very coarse radar resolutions of 1 km by 1 km and 
2.6 km by 3.1 km gave the best results for wet moisture conditions 
while a finer resolution of 93 m by 100 m was found to yield superior 
results for dry to moist soil conditions. 
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ABSTRACT 


Simulation techniques were onployed to generate orbital radar 
imagery of a 17. 7"km by 19«3“km test site near Lawrence, Kansas. The 
simulations were produced for a radar operating at ^.15 GHz, with HH- 
polari zati on , and over an image swathdef ined by an angle of incidence 
between 7.5“ and 9.3“ from an orbital altitude of 600 km above mean 
sea level. Images were simulated for synthetic aperture radars (SARs) 
with resolutions of 20 m by 20 m, 93 m by 100 m, and 1 km by 1 km, and 
also for a real -aperture radar (RAR) with a resolution of 2.6 km by 3.1 km. 
For each simulation, the power received by the radar was determined with 
empirically derived formulae from a static terrain model and a dynamic 
meteorological and agricultural model that acted on each of the 20-m by 
20-m pixel elements within the test site. The terrain model included 
surface elevation, land-cover category, and A-horizrn soil texture. The 
dynamic model established dally distributions of surface- layer soil moisture 
from rain gauge and pan evaporation data using a water-budget approach 
dependent upon canopy cover, stage of crop development, soil type, and 
a Monte Carlo synthesis of between- and within-field variability. Radar 
simulations were performed at five-day intervals over a tventy-day period 
and used to estimate soil moisture from a generalized algorithm requiring 
only received power and the mean elevation of the test site. The results 
demonstrate that the soil moisture of about 90 per;c.nt of the 20-m by 20-m 
pixel elements can be predicted with an accuracy of ± 20 percent of field 
capacity within relatively flat agricultural portions of the test site. 

Radar resolutions of 93 m by 100 m wi tn 23 looks or coarser gave the 
best results, largely because of the effects of signal fading. For 
the distributions of land-cover categories, soils, and elevation in the 
test site, very coarse radar resolutions of 1 km by 1 km and 2.6 km by 
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3.1 km gave the best results for wet moisture conditions while a finer 


resolution of 
to moi St soi 1 


93 m by 100 m was found to yield superior results for dry 
condi tions. 
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t.O INTRQDIJCTIQN 


(mag« simulation tochniquas offer an Ideal mecnanlsm by which to 
examine radar resolution requiremenis for characterization of area 
extensive terrain features* Radar Images have been real 1st I cal ly 
simulated for the Seasat L^band radar and aircraft Imaging radar 
systems when given appropriate terrain models and statistically 
accurate radar cross-sections for constituent scene elements [26,27], 

The validity of the Image simulation approach In defining optimal 
sensor resolution for a given application rests. In large part, upon 
the realism of the terrain model with respect to the application. 
Since our concern rests with an accurate assessment* of near-surface 
soil moisture, special emphasis must be placed on Incorporation of 
credible Instantaneous soil moisture distributions at a scale less than 
or equal to that of the sensor resolution and the modeled moisture must 
be distributed over relatively large extents (If the entire radar Image 
swath Is to be simulated). This Is certainly a non-trivlal objective. 

The wel l-documen-** j variability of field measured soil properties 
[17,18,19.20], especially In the dynamic surface horizon, even within 
"homogeneous” test plots precludes the use of area-extensive point 
measurements of surface soli moisture as adequate Inputs Into a data 
base for radar image simulation. Thus, it seems reasonable to use a 
modeling approach to define the surface- layer soli moisture of a given 
terrain element as a function of time, with the model based upon the 
physical properties of that element. 

This repv'rt will briefly review the results of earlier satellite 
radar simulation studies that attempted to define resolution 
requirements In soil moisture estimation and will Interpret these 



findings with respect to the procedure used to model soil moisture. 
Section 2 will develop a more detailed and realistic soil-water model 
which will feed a series of satellite radar simuiatio..s presented and 
analyzed in Section 3. These simulations will Include three synthetic 
aperture radars (SAR) with nominal resolutions of 20 m, 100 m, and 1 km 
and also a real -aperture radar (RAR) with a nominal resolution of about 
3 km. 


1.1 Erior Badar 5iiaulatlafl.fifiSMi.s,fQr..Aas.as5infiRt. Qf 

An earlier examination of resolution requirements of an orbital 
SAR used for estimating surface soil moisture by image simulation 
techniques Cl >23 made use of the radar configurations shown In Table 

1.1. The satellite, operating at a frequency of 4.75 GHz with HH 
polarization. Is given an orbital altitude of 600 km above mean sea 
level and an antenna pointing angle range from 6.4 degrees to 20.0 
degrees which results In a nominal ang I e-of- I nc I dence range of 7 
degrees to 22 degrees, assuming a curved-earth surface as shown In 
Figure 1.1. Of the fourteen system design options presented In Table 

1.2, three synthetic-aperture radar* were chosen for simulation - Cases 
2, 4, and 8 - which yielded approximate resolutions of 20 m with 12 
looks. 100 m with 23 looks, and 1 km with 1,000 looks respectively. 

The radar simulations operated on a data base constructed at a 
20-m X 20-m grid cell resolution for a predominantly agricultural 
region (17.7 km x 19.3 km) located just southeast of Lawrence, Kansas 
Cl»2j. The net percent area occupied by each target class Is shown In 
Table 1.3. In addition to target class as Interpreted from U-2 color 
IR photography, each of the 800,000 grid cells In the data base was 
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TABLE 1.1 


Spacecraft Radar Configuration 


Spacecraft Height 

h 

600 km 

Radar Frequency 

f 

4.75 GHz 

Angle of Incidence Range 

0 

70 . 22 “ 

Antenna Pointing Angle Range 

6 

6.39“ - 20.02° 
(curved earth) 

Receiver Noise Figure 

F 

4 (6 dB) 

Minimum Signal-to-Noise Ratio 


4 (6 dB) 

System Loss Allowance 


2 (3 dB) 

Minimum Scattering Coefficient 

0 

a • 
mm 

- 21 dB 

Antenna Length 

D 

8.7 m (8.7 m and 15 m 
for RAR and 5 m for 
RAOISAR) 

Antenna Radiation Efficiency 

n 

0.75 

Ground Swath-Width 


143 km 

Slant Swath-Width 


38.5 km 









Potential System Design Options 
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TABLE 1.3 

Area Percent of Total Data Base 
As s iqned to Each Target Class 


Target Class 

Percent of 
Total Area 

Roads 

3.76 


Railroads 

0.12 


River Bridges 

0.01 


City Structures 

0.85 


Rivers 

2.19 

2.67 

Lakes, Ponds, 
Impondments 

0.48 

Smooth Bare Soil 
(RMS height < 2 cm) 

6.63 


Medium Rough Bare 
Soil 

(2 cm < RMS height 
4 cm) 

4.92 


Rough Bare Soil 
(RMS height > 4 cm) 

2.75 


Mown Pasture 

7.06 


Pasture 

15.93 


Alfalfa 

4.15 


Wheat 

6.65 


Sandbars 

0.35 


Deciduous Trees 

13.03 


Soybeans N/S Rows 

5.46 


Soybeans E/W Rows 

5.62 

11.08 

Milo N/S Rows 

2.76 


Milo E/W Rows 

2.27 

5.03 

Corn N/S Rows 

8.32 


Corn E/W Rows 

6.61 

14.93 


* N/S refers to crops planted with rows running north to south. 

* E/W refers to crops planted with rows running east to west. 
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assigned an elevation from USGS 7.5-mlnute-serIes topographic maps and 

I 

a surface soil texture from USOA/SCS county soil surveys. 

A set of four general soil moisture conditions was computed for 
this data base for subsequent Integration Into the target^sensor 
Interaction submodel of the radar simulation program. Surface soil 
moisture was calculated by a highly simplistic dally accounting of 
Incident rainfall and evaporative losses with some limiting assumptions 
so that soil moisture could not exceed saturation or be less than the 
hygroscopic coefficient. A storm submodel generated precipitation with 
a very broad Gaussian distribution of Intensity across a linear storm 
track. The water budget model did not account for the fol lowing: 
•surface slope 

•effects of soil type on soil bulk density or hydraulic 
conductivity 

•transpiration by a vegetation canopy 
•soil profile effects 

•dependence of hydraulic conductivity on soli moisture 
•presence of standing water on surface. 

A summary of the postulated soli moisture conditions Is presented 
In Table 1.4. The soil water-budget model produced 0-5-cm soil 
moisture values as volumetric moisture (g/cm^) which were subsequently 
converted to a percent of field capacity (MpQ) basis by an empirical 
expression dependent upon soil textural components. 

Elevation Information was used to calculate the local angle of 
Incidence 9^^ from range and slope Information for each pixel. The mean 
backscattering coefficient for a cell was computed from a series of 
empirical relationships found to be statistically valid for each 
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Values of Rainfall and Evaporation Used to Derive the Hypothetical Soil Moisture Conditions 
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combination of Mpc, and target category CO. For a given sensor 

-0 

configuration and terrain element area» c of a pixel was converted 
Into an averaged received power Pp by the radar equation. Since, for 
the three resolutions simulated, the number of looks N was always less 
than Infinity, Pp was randomized for Rayleigh fading with the 
appropriate number of looks to a received power Pp at the antenna. 

The resultant Images of Pp were then machine processed Into soil 
moisture estimate maps using only celt position and mean elevation of 
the whole 17.7-km x 19.3‘-km data base (from which an estimate of ^nd 
area could be deduced for each pixel assuming a spherical earth). The 
moisture content of each cell was then estimated from 

Mpc = - f(0)3/g(e) (1,1) 


where 


MpQ = estimated percent of field capacity of the 0-5-cm layer, 
a° = scattering coefficient (In dB) as estimated from Pp 
and range (assuming spherical earth), 

9 » angle of Incidence estimated from range to target, and 

f(6) and g(e) = empirically determined polynomial expressions 

[ 1 ]. 


Finally, the radar soil moisture estimate accuracy of a given 
simulation was deduced by comparing the map of MpQ to the spatial 
distribution of Mpc defined by the soil water-budget model . This 
process yielded an Imperfect comparison since the radar Imaging process 
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significantly distorted the geometry, and positioning errors for even 
small hills in the data base often lead to offsets of 100 m between 
and where i and J are relative positioning coordinates. 

This entire process and the data and assumptions on which it is based 
are treated in a more comprehensive fashion by Ulaby, et ai. Cl end 
23. 

Results of the above comparison procedure for each of the 20-‘m x 
20-m grid cells In the actual moisture distribution generated by the 
water-budget model are shown in Figures 1.2 and 1.3 for the whole 
800,000-cel I data base and the 183,000-ceil river floodplain region 
respectively, for moderately dry to moist soil conditions. Results 
such as these, between angles of 7 degree and 20 degrees, showed that 
estimates of soil moisture from the 100-m and 1-km resolution radar 
systems yielded significantly less estimate error than that afforded by 
the 20-m resolution radar. Furthermore, if nonagri cultural categories 
such as buildings, roads, bridges, and water (where soil moisture Is 
undefined on the map) are excluded from the comparison, and the 
comparison is also limited to the relatively flat floodplain of the 
river (where positioning errors are minimized), between 81]t and 90K of 
the 20-m x 20-m grid ceils are found to have less than 20 % absolute 
error in the estimate in Figure 1,3. it shouid be noted here that 
these comparisons are performed on a 20-m x 20-m basis, hence ail 

sensor resolutions are being evaluated as though they are predicting 

moisture at a 20-m scale. 

The effect of averaging the Input moisture distribution Mj.^. to a 

t-km x 1-km scale is shown in Figure 1.4. Also shown are the 

distributions of predicted moisture resulting from Interpretation 
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Figure 1.2. Soil moisture estimate accuracy ^or the entire simulation 

data base with moisture condition 3, ten days after thunder- 
storm, and at anoles of incidence between 7.5° and 9.3°. 
Absolute difference is comouted at a 20 meter resolution from 

r2]. 
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Figure 1.3. Soil moisture estimate accuracy for agricultural categories 
on the floodplain with moisture condition 3, ten days after 
thunderstorm, and at angles of incidence between 7.5* and 
9.3“. Absolute difference is computed at a 20 meter resolution 
from [2]. 
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Moisture Condition 3: 10 Days after Thunderstorm 
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30 40 50 60 70 80 90 100 110 

Percent of Field Capacity 


Comparison of actual soil moisture distribution {as 20 m x 20 m 
and 1 km X 1 km averages) with the distributions predicted from 
simulated radar imagery at 20 m and 1 km sensor resolutions. 


of the 20-m resolution SAR and the 1-km resolution SAR. Obviously, the 

A 

from the I -Km SAR more close I tracks the actual distributions of 

over most of the range. The combined effect of varying both radar 

resolution an.^ the scale of comparison (or resolution of the moisture 
estimate) Is demonstrated by Figure 1.5. This figure rescales the 
results plotted on Figure 1.2 for SAR resolutions of 20 m and 1 km (as 
used to estimate soil moisture at a 20**m resolution) and shows that an 
additional "Improvement" In accuracy Is achieved when the i km SAR Is 
used to estimate moisture at a 1-km x 1-km scale. This apparent 

Improvement Is not unexpected and Is merely a consequence of the 

reduction of local variance In average Mp^. 

A summary of results for al I four moisture conditions generated 
from the water-budget model Is shown In Figure 1.6. Here the percent 

of 20-m X 20-m grid cel Is with an absolute moisture estimate error of 

less than 30)( of field capacity Is plotted as a function of resolution. 
The trend toward Increasing estimate accuracy by using coarser 

resolution sensors Is apparent for all moisture conditions a 5 thouph 

this trend seems to break down for extremely dry soil moisture 
conditions. 

These findings have been attributed to several considerations C23» 

DM as generated by the water-budget model varies dramatically 
F C 

only between adjacent soil types (typically a soil type Is on the order 
of 100 m to 1 km wide); 

2) the coarser resolutions act as low-pass spatial filters and 
average the local "noise" effects of slope, canopy cover type, row 
direction, and surface roughness, and 

3) the effects of water bodies, cultural targets, and forested 
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Resolution (A/leters) Resolution (Meters) 

FigurG 1.6. PGrcent of data base where estimateo moisture Mpr is within ±30% of 
actual moisture Mpc. Angle of incidence 0 is 7.5° to 9.3° relative 
to nadir. Accuracy level 3 is approximately ±30% of field capacity. 


areas are averaged over much larger regions by a coarse resolution 
sensor. 

Implicit in these conclusions la that the most appropriate 
(accurate) sensor resolution for moisture mapping Is Inextricably 
linked to the scale of variance In target parameter of Interest: 
soil moisture. ThP: Is further complicated by the density of the 
distributions of observable "noise" parameters such as cultural 
targets, water bodies, hills, and within-fleld surface roughness 
effects. Hence, It Is possible that a different set of assumptions 
regarding the nature of actual moisture distribution may lead to very 
different conclusions. 

Figure 1.7 shows the cumulative distributions of soil moisture 
MpC produced by the simple water-budget approach described above. 
Whi I e the net moisture conditions are very different, ranging from 
saturated soil conditions on Day 5 to near The hygroscopic coefficient 
or Da” 35, only Day IS (Moisture Condition 3) encompasses a wide range 
of Instantaneous soil moisture conditions. For Day 15, the 
discontinuities present on the cumulative graph demonstrate that 
variance In moisture Is primarily controlled by soli type and only 
weakly by the minor variance In rainfall from the storm submodel. 

1 .2 Modificatio ns of the Simulation Data Base 

As a consequence of the simplicity of moisture distributions 
produced by the above water-budget model, the conclusions reached Cl »23 
are reevaluated using a more complicated and realistic accounting 
procedure as well as Introducing sources of between- and within-field 
variance other than cover category (Table 1.3) and soil type. 
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Figure 1,7, Cumulative frequency distributions of "actual" soil 
moisture for the four simulated conditions generated 
from the simplified environmental model given in Table 1,3, 
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respectively. Figure 1,8 presents a diagram of the simulation process 
modified by Section 2 and used to generate the Imeges analyzed In 
Section 3. 

Significant Improvements In the data base as a dynamic model for 
radar simulation studies of agricultural terrain Include: 

•Incorporation of cropping calendars to al I om for 
additional between- field but with In-crop variability; 

•addition of within-fleld variability In moisture 
by a Monte Carlo approach; 

•a more comprehensive soil wattr accounting model (SWAM). 

A diagram of surface conditions relevant to SWAM is shown In 
Figure 1.9. The model basically consists of three potential layers: 
standing water, the upper 5 cm of the soil, and a crop root zone 
(nominally extending to a 1-meter depth). The water table is always 
treated as being well below 1 meter. The model, while still 
simplistic, treats both infiltration and evapotranspi ration as 
dependent on dynamic soil and crop canopy conditions. In addition, a 
more robust storm model leads to a greater variability in incident 
rainfall, which Is subsequently allowed to run off or to Infiltrate 
depending upon the intensity of the rain, antecedent soil moisture, and 
local surface slope. 
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Figure 1.9. Local environmental conditions pertinent to the soil water accounting model. 



2.0 DYNAMIC SOIL WATER ACCOUNTING MODEL (SWAM) AND VARIANCE 

iN,SfllL, MOISTURE 

The purpose of a soli water-budget mode! within the context of 
real 1st I c radar image simulation is to generate a distribution of 
near-surface (0-5 cm) soil moisture conditions at the spatial scale of 
the static terrain data base (20 m x 20 m) which responds to both 
static conditions (soil type, cover type, and surface slope) and 
dynamic conditions (crop stage, rain, and potential evaporation) on a 
time scale relevant to both the dynamics of the process and the orbital 
mechanics of an Imaging satellite (dally basis). While many excellent 
water-budget models are available for various applications In agronomy 
and hydrology C3 to 7j, no single model meets all the above criteria. 
Indeed, most such models require more detailed Information on soli 
profile characteristics and weather conditions than Is readily 
available for the simulation area and are designed to operate at a 
level much less than field size over time Increments significantly less 
than one day, or conversely, they are most appropriately applied to 
very coarse Integration times on the order of weeks for a simple sot of 
Input parameters and at a macroscopic level much larger than field 
size. 

Because of the large size of the data base (approximately 800,000 
grid ceils) it Is necessary to tailor a model that emphasizes the 
surface horizon, requires a minimum of Information as to soil profile 
and detailed local weather conJItlons, and yet Is still sensitive to 
dally variation In soil moisture. A schematic of the final process 
model Is shown In Figure 2.1; It consists largely of the following 
components : 
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Figure 2.1. Dynamic Soil Water Accounting Model (SWAM). 
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•storm model , 

•surface runoff model, 

•crop development submodel, 

•evapotranspiratlon model, and 
•an Interlayer redistribution model . 

When given dynamic Inputs of crop type, crop stage of development, 
rainfall, and potential evaporation, the model acts upon the static 
terrain model to yield dally projections of 0-5-cm soil moisture for 
each grid cell. It also governs the redefinition of canopy cover 
categories based on crop-calendar changes or local flooding conditions, 
and these categories are then used as Input to the radar Image 
simulation program's target-sensor Interaction model. 

2.1 Storm Model 

The storm model used In SWAM combines a rain Intensity 
distribution function determined by storm type with rain-gauge data and 
a Monte Carlo appraoch to storm-track positioning. For computational 
ease, the simplifying assumptions are made that all storm cells proceed 
along an east- west axis across the data-base matrix and do not begin or 
terminate (temporally) within the data base. 

Dally rainfall In cm as recorded at Clinton Reservoir, Kansas In 
1978 Is shown In Figure 2.2 for a 220-day period as extracted from NOAA 
monthly summaries of local climatological data. For each rainfall 
event, the duration of the rainfall was also recorded and average 
rainfall Intensity was computed from 
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Reservoir, Kansas during 1978. 


ig • Rg/t 


( 2 . 1 ) 


where 

Ig ■ rainfall Inteiiilty at the gauge, cm/hr 
R • net dally rainfall, cm 
t ■ storm duration, hours. 

Each storm was also classified as a local convective event, as part of 
a largo frontal event, or as Intermediate. This three-scale storm 
classification serves as the basis for rain Intensity functions that 
are based, in turn, upon results which Indicate that Intensity decays 
exponentially from the peak Intensity to some minimum value at the 
perimeter of the ra In-producing storm cell C9-123. The function takes 
the form: 

1(D) » l„ax * [At + By * expCCy * Y)]/D (2.2) 

where l^^^ » Intensity at the storm center, 

A,B,C * are fitted values for each storm type T 
(based on C)0,1 

Y ■ the lateral radius of the storm cell, and 
0 ■ distance between I and I max* 

The range of fitted values from convective to frontal -type storms are: 
10.8 to i»5.7 for A, 

163.5 to 90*8 fof and 

-0.20 to -0.09 for C. 
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since the center of the storm cell Is not known relative to the 
measured gauge data. It Is synthesized by an equl<>probabl 1 Ity 
random-number generator. Intensity along the center of the storm 
trajectory Is then computed by Inverting Eq. 2.2 to the form 

•max ■ Ig * D/[Ay + By * exp(Cj * Y)], (2.3) 

where D ■ random distance between the storm track and the rain gauge. 

A second randomization Is then used to position the storm track 
with respect to the 800,000 grid cells In the data base. As a result, 
it Is possible for specific storm events measured on the gauge data to 
completely miss the simulation area; however, the application of 
appropriate limits to the random-number generators makes this a rare 
occurrence. 

Thus, given the maximum rainfall Intensity, the trajectory of the 
storm cell, the storm type, and the extent or the storm, the rainfall 


Intensity of any grid cel 

1 can 

be calculated from Eq. 

(2. 

2). For 

the 

20-day portion 

of the 

rain 

gauge data after Jul 

Ian 

day 153, 

the 

Intensity of 

al 1 storm 

events Is plotted as 

a 

function 

of 


north-to-south distance within the simulation data base In Figure 2.3. 
A distance of zero designates the northern edge of the data base. 
Three rain events occur within this period In June, during which the 
simulated radar overpasses occur on a five-day revisit period. Two of 
the storms are local convective showers, both having measured rainfall 
of 0.8 cm over one-hour durations. The randomization of storm-track 
location separates the two events In space and gives the one on Jul lar 
day 156 a slightly higher maximum Intensity. The rainfall on Julian 
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Figure 2.3. Rain intensity over the simulation area q« modeled from rain-gauge data 
during the 20-day period covered by the radar simulations. 


day 172 Is much mora extansiva and has a duration of alg^t hours. Tha 
nat dally rainfall as distrlbutad by thasa thraa ovants ovar and around 
tha simulation araa Is shown In FIgura 2.4, Tha maximum Incident 
rainfall ovar tha 20-day parlod Is saan to ba approximately 9 cm, which 
Is raasonabla for tha simulation rag Ion In mId-Juna. 

2.2 Surfaca Runoff Modal 

Tha s*j»-faca runoff modal considars only tha nat affact of local 
surfaca slopa and doas not axpticitly account for water retention and 
Impoundment by soli surfaca roughness, tillage practices, and tha 
presence of terraces. The water available for drainage as lateral 
surfaca flow Is equal to tha sum of standing water remaining from tha 
previous dally accounting parlod plus tha Incident rainfall In excess 
of that which can Infiltrate the surface layer of the soil during tha 
rainfall event. Tha Infiltration ra+o Is linked to the antecedent 
moisture of the surface layer. Thus the potential drainage Op becomes i 

Dp - SW + t«(l - K) (2.4) 


where 

SW • standing water 
t ■ duration of rain event 
I • IntansI ty of rainfal i 

K ■ hydraulic conductivity of the surface layer. 


The actual drainage 0^ I:. computed from Dp and local surface slope 


Net Rainfall, cm 
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Figure 2.4. Net daily rainfall over the simulated area during the 20 day 
period covered by the radar simulations. 
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Da = Dp » (1,0 - 0.8“) 


(2.5) 


where a* the slope angle of the surface from horizontal In degrees. 

The term 1,0 - 0,8°‘ is defined as the drainage coefficient and Is 
plotted versus surface slope (In percent) In Figure 2,5. While most of 
the slopes in the data base are less than 6 percent for agricultural 
fields, local slope can exceed 15 percent for some of the hills 
adjacent to the Kansas River floodplain. 

2.3 Evapotransplratlon Model 

Evapotranspiratlon Is calculated differently for cropped and bare 
soil surfaces. For bare soil surfaces the actual evaporation Is 
depleted solely from the soil surface layer, while for vegetated 
surfaces a static root distribution model removes 20 percent of the 
actual evapotransplratlon from the 0-5-cm layer and removes the 
rema'ning 80 percent of actual evapotranspi ration from the "root zone." 
For simplicity, the "root zone" Is assumed to be one meter In depth and 
Is treated as a constant with time and for all crops. 

For bare soil, actual evaporation, AE, Is computed from potential 
evaporation, PE, as limited by antecedent soil moisture In the surface 
layer and soil hydraulic properties. Accounting Is performed on a 
dally basis using the daily pan evaporation recorded at Clinton 
Reservoir, Kansas and Is shown In Figure 2.6 for 1978. 

An experimental nrodel Is used to calculate actual evaporation from 
potential evaporation PE Cl3j: 
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Figure 2.5. Variation in drainage coefficient as a function of soil surface slope 
in the soil water accounting model (SWAM). 
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rvoir, Kansas during 1978, 


AE » PECA + B(MR) + C(MR)2 + 0(MR)33 


( 2 . 6 ) 


where A, B, C, and D are empirically derived coefficients dependent 
upon PE Cl33» and the moisture ratio MR Is given by 

MR = (SM - WP)/(FC - WP) (2.7) 

where SM » the measured soil moisture, 

WP = the wilting point of the soil, and 
FC = the field capacity of the soil, and 0 ^ MR £ 1.0, 
and the potential evaporation PE Is 

PF = kp » Epan» 

where kp = pan coefficient and 

Epan = measured pan evaporation. 

Regression fits of experimental data yield Cl33: 


A » -0.05 + 0.732/PE (2.8) 

B= 4.97 - 0.661 PE (2.9) 

C = -8.57 t 1 .56 PE (2.10) 

0 » 4.35 - 0.88 PE (2.11) 


Thus, for a given day, all terms In Eq. 2.6 are constant for all grid 
cells except the moisture ratio, which Is dependent on the antecedent 
soil moisture and the gross water-retention characteristics of each 
soil. 

For vegetated soils, the actual evapotranspi ration, ET^pop, is 

3 ^* 



computed by a modification of the Blaney-Crlddle formulation used In 
estimating crop Irrigation requirements [8,20. Although the method Is 
designed for effective Integration periods of weeks to months, the 
simplicity of Its Input requirements makes this a practical approach 
for such a large number of grid cells. Basically, crop consumption of 
water over the rooting depth varies with temperature, length of day, 
available soil moisture, crop type, crop stage of growth, relative 
humidity, and windspeed. To simplify the formulation, average measured 
values of temperature, day length, relative humidity, and windspeed are 
assumed on a seasonal basis for the simulation area. The resultant 
expression for ETcpop becomes [83s 

ETcrop = kc * ETq (2.12) 

where kc = crop coefficient 

ETo = potential evapotranspiratlon. 

While ETq Is generally computed on a mean monthly basis by the Penman 
equation, Eq. (2.12) becomes more sensitive to <^a!ly changes In 
root-zone soil moisture fo^* ETq = AE as computed by Eq. (2.6) for the 
rooting zone. 

Crop coefficient as adjusted for mean local climate Is plotted In 
Figure 2.7 as a function of the number of days after planting for 
several of the crop covers Included In the data base. Crop consumption 
of water Is seen to be dependent on both crop type and stage of crop 
development. 
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Figure 2.7. Crop transpiration coefficient as a function of time from planting for several 
representative crops found within the simulated area. 



2.4 Crop PeyglQPinent Model 


The lengths of time required for a given agricultural field In the 
simulation data base to progress from one crop-development stage to the 
next Is established from data gathered by the Statistical Reporting 
Service of the United States Department of Agriculture, Figure 2.8 
shows the crop reporting districts. The simulation area lies at the 
northern limits of the East Central reporting district of Kansas (No. 
6). Figure 2.9 presents a summary of mean crop development over a 
1‘0-year period as enumerated by AgRISTARS Cl 53 for this crop reporting 
district. These percentages are used to define crop development stage 
within the simulation on a f leld-by-f leld basis. Thus, each distinct 
agricultural field In the data base Is assigned one of ten codes which 
allows that field to be Individually allocated to one of 10 planting 
dates. Hence, there are ten different absolute crop calendars possible 
for each crop type Identified In Table 1,3. 

Planting dates are randomly assigned to field codes for a specific 
crop based upon Figure 2.9. This procedure results In the Introduction 
of a significant source of between- f I el d variance of soil moisture 
within a given crop type due to the effect of crop development stage on 
evapotranspiratlon. It also allows for a given field to have Its 
target classification changed In Table 1.3, since a medium-rough bare 
field becomes a cropped field after emergence, and finally reverts to 
bare soil status after harvest. As Implemented, this procedure gives 
tiie data base a dynamic crop-category mix that can be modified to match 
regional agricultural practices such as double cropping or dynamic soil 
surface roughness conditions. 
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CROP-REPORTING DISTRICTS 
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2.5 


Infiltration of water Into the surface layer, percolation of water 
into the root zone, and capillary recharge of surface layer moisture 
are controlled by the water-potential profile as limited by soil 
structure. Applying Darcy’s formulation to the flow In unsaturated 
sol Is yields: 

V - K(9) * VH (2.13) 

where 

V ■ the flow velocity, 

K(9) » the hydraulic conductivity as a function of volumetric 
soi 1 moisture 9 , and 

7H * the gradient of the hydraulic head H. 

For the soUc Included In the data base, the volumetric moisture 9 Is 
approximately known only at matric potentials of 1/3 bar and 15 bars, 
and no ready data source exists for K(6) except at saturation. 
Therefore, the approach given by Eq. (2.13) for laboratory 

measurements or modifications based on field measurements C63 Is not 
tractable. However, a rough first-order approximation of the effect of 
K(6) on flow rate V Is Incorporated which allows K to decrease in a 
linear fashion from that measured at saturation as given by; 

9 “9 

V - Ksat • (1.0 - ) (2.14) 

“max 

where Kggt “ measured saturated hydraulic conductivity 
®max “ porosity, and 
"max * ' " 

Ao 


Pjj ■ soil bulk density, g/cnP 
Pfj ■ soil specific density, g/cm3 

Assuming that a! 1 the soils in the data base have a specific density of 
2.65 g/cm^, Eq. (2.14) reduces to: 

V « 2.659Ksat/(2.65 - P^) • (2.15) 

It is recognized that this approximation can seriously overestimate 
K(9) of clays at low moisture contents as shown In Figure 2.10. 
However, since the function is uniformly applied to Infiltration at the 
surface and to percolation Into the root layer, the net effect of this 
error on surface-layer moisture will be reduced af te ' several 
accounting periods Inasmuch as water will readily drain, and It is the 
surface water content that drives the radar simulation model. 

In addition, several limiting conditions are placed on percolation 
and capillary recharge. Percolation only takes place when 
sun ace- layer moisture exceeds field capacity and only to the extent 
that surface- layer moisture Is reduced to field capacity. Also, 
capillary recharge is only considered as significant when 
evapotransplratlve losses have reduced the surface-layer soil moisture 
below the wilting point, thus allowing for the formation of a 
discontinuity In moisture at the layer boundary that Is related to the 
formation of a surface crust. 

Whi'e the above approach Is not rigorous. It allows for the ready 
calculation of approximate flows with the proper directional Ity ano 
requires minimal Inputs and computation for the very large number of 


U] 


ORIGINAL cj 

OF POOR QUALITY 



Figure 2,10. Hydraulic conductivity as a function of volinetric 
soil moisture. 
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grid celts In the data base. An alternate approach can be used 
for soils where K(e) Is measured for each soil at a low value of 
moisture: 

K(0) = Ksat * expCa » Ce-9maxJ] (2.16) 
where a Is a free parameter set to fit the measured values. 

2.6 Within-Field Variab ility In Surface Soil Moisture 

Prior to radar Image simulation, the surface layer soil moisture 
values determined by SWAM for each 20-m x 20-m grid cell are randomized 
to approximate the natural variability In soil moisture measured within 
'•homogeneous” fields. Rather than rely on field-based studies of 
moisture variance for various plot sizes as computed from 
point-sampling techniques C17-20], a variance estimate based on areal 
measurements of soil mcisture at a scale close to 20 m x 20 m is deemed 
more appropriate. 

The areal surface-moisture sampling Instrument considered Is the 
4.75-GHz airborne scatterometer flown by NASA/JSC over a test site near 
Colby, Kansas In July and August 1978. The test site consisted of 
forty-two 40-acre agricultural fields selected because of their 
expected wlthln-fleM homogeneity. While the 260-km2 test site was 
overflown seven times, data from only the first six flights are used In 
this analysis; the last flight produced acceptable scatterometer data 
but suffered from a camera malfunction which makes ground referencing 
of the backscatter data uncertain. The 4.75-GHz fan-beam scatterometer 
was flown at a mean altitude of 460 m ana the data was processed by 
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NASA/ JSC and averaged over 0.5-second time Intervals. At a typical 
aircraft ground speed of 278 km/ hour, this yields a reported 
scatterometer fore-aft ground resolution of 37 m. Since the 
cross- track 3-dB beamwidth Is 3 degrees for the 4.73-GHz scatterometer, 
at a 10-degree angle of Incidence, the cross- track resolution becomes 
25 m. The radar backsca'ter Ing coefficient a° from this 37 -m x 25-m 
resolution element Is the mean of 90 Independent samples from frequency 
averaging as calculated by 

Nf = 2(AD)2 cos3e/Xh (2.17) 


where: 

= number of independent samples produced by frequency 
averaging, 

AD * fore-aft ground resolution, 

9 « angle of incidence, 

h = aircraft altitude, and 
\ = wavelength. 


For 90 independent samples, the uncertainty in a° due to fading Is 
less than +/- 0.15 dB assuming Ray(elgh fading statistics; hence, the 
scatterometer* s measurement precision can be considered to be quite 
good. 

For the 40-acre test fields at Colby, the variance In 0-5-cm soli 
moisture has a mean value of 6.5 percent on a gravimetric basis from 45 
sampling locations In each of 154 field observations [223. Ignoring 
fields with crop row direction orthogonal to the aircraft Might path, 


the least-squares linear correlation between 0-5-cm soil moisture and 
a°(dB) Is found to be 0,80 when comparing field mean values of with 
moisture at an Incidence angle of 10 degrees and HH polarization 
and 0.92 (FIgurj 2.11) when comparing on a subfield basis with the 
nearest point measurements of soil moisture C23] at an Incidence angle 
of 15 degrees. While both of the above sets of linear correlation 
results are based upon only the first two flights, the addition of data 
from all six flights does not significantly alter the strength of the 
linear correlation between a° (dB) and near-surface moisture. In a 
comparison of f’eld averages of moisture and a° , the linear correlation 
p Is found to vary between 0.82 at 10 degrees and 0.87 at 15 degrees. 
Due to the variability In sampled soil moisture found even within 
small plots, and due to sampling error Cl7 to 20l]» the maximum linear 
correlation coefficient attainable Is limited, even for a perfect 
moisture sensor, by the variance In point-sampling error. For 
point-sampling error distributed with a zero mean and a variance of 
, and assuming moisture and sampling error to be Independent, It can 
be shown that C243 


P (p 


"SM> • 


1 + 


n 


( 2 . 18 ) 


where 

M = mean soil moisture estimate, 

SM 

0^ = sampling error variance, and 

n 

^cu ~ true soil moisture population variance. 
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ALL FIELDS, EXCEPT WITH PERPENDICULAR ROWS, WERE USED 

R-SQUARE «• 0.846 


Equation (2.18) shows the correlation coefficient to be Inversely 
proportional to sampling error variance and directly proportional to 
the true variance of soli moisture. Thus, assuming that the true soli 
moisture variance Is that measured within the Colby test fields, = 
6,5%, and assuming that the sampling error variance Is within the range 
found for dry soil sample weights less than 100 grams C203, 
i.2,0%, then the maximum linear correlation coefficient calculated from 
Eq. 2.18 Is 0.87^, P S.0.94. As a result. It Is not statistical ly 
unreasonable to assume that the scatterometer Is a perfect moisture 
mapper. 

Assuming for the moment that the Intrafleld variance In crop 
canopy conditions and random surface roughness Is smal I enough to be 
negligible, and assuming that the variance In radar backscatter at 4.75 
GHz, HH polarization at a 10-degree angle of incidence Is solely 
dependent in a linear fashion on surface soil moisture, the Colby 
4.75-GHz aircraft scatterometer data can be used to define an estimate 
of the natural variance In the true soli moisture population between 
sensor resolution cells of 37 m x 25 m. Figure 2.12 shows a measure of 
wlthln-fleld variance In CY, plotted as a function of a°. 

SD 0 

CV ^ 10 log (1 + —^) (2.19) 

- 0 
a 

variance coefficient, 

within-field standard deviation in , m/n, and 
field mean o" , m/m. 


where; 

CV 

SD 0 = 

a 

-0 
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Figure 2.12. Variance in aircraft radar backscatter due to intrafield variability of moisture 

for "homogeneous" test fields at Colby, Kansas in 1978. Scatteroneter characteristics 
are 4.75 GHz frequencies, HH polarization, 10° angle of incidence, and a nominal 
resolution of approximately 30 meters. 


For sample sizes between 6 and 10» the within field variance of a° 
Is seen to be Independent of which, for a perfect soli moisture 
mapper. Implies that CV is independent of soil moisture (which ranged 
from the wilting point to saturation over the course of the six flights 
at the Colby site). The mean value of CV Is *ound to be 0,78, which 
corresponds to an average within-fleld uncertainty In of + 0.70 dB 
to - 0.95 dB. Application of these uncertainties to the empirical 
smooth bare soil regression equation between (dB) and 0~5*cm soil 
moisture 


<^°(dB) = 0.149 Mpc - 15.49 


( 2 . 20 ) 


where: 

Mpc “ 0-5 cm percant of field capacity. 


yields an uncertainty In Mp^ of + 5.3 to - 6.4)1. 

The output of the SWAM budget model for the 0-5-cm soli layer Is 
randomized on a grid-cell basis by a Gaussian random-number generator 
with a standard deviation of 6 percent Mp^. This within-fleld moisture 
variance Is added by pixel for each moisture distribution prior to 
calculation of the mean grid cel I a° in the target^sensor interaction 
model for radar Image formation. 

An Image format example of the Impact of the SWAM plus 
within-fleld variance In 0-5-cm soil moisture distribution Is shown In 
Figure 2.13. The figure compares Moisture Condition 3 from the simple 
water budget CO with the moisture distribution present over the same 
3-km X 5-km portion of the data baso as produced by SWAM on Julian day 

ks 
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Moisture Condition 3 from the simplified environmental model 


(b) Julian day 1 63 from the dynamic soil water accounting model (SWAM). 


Cotnparison of the spatial distributions of soil moisture 
resulting from a) the simoMfied environmental model 
given by Table 1.3, and b) the dynamic soil water accounting 
model. The general soil moisture conditions are similar 
for both models and grey tone is proportional to moisture. 






TABLE 2.1 


Comparison of Mean Soil Moistures from the Simplified Environmental 
Model with those from the Dynamic Soil Water-Account Ing Model 
for the Same Agricultural Fields Under Similar 
Regionpl Moisture Conditions 


Fi el d^ 

Soil 

Crop 

Moi sture 
Map V 

Mean % 
of Field 
Capaci ty 

■ 

Number 

of 

Pixels 

A1 

Loam 

Pasture 

MFC 3 

32.0 


36 

A1 

Loam 

Pasture 

Day 168 

97.17 


36 

A2 

Silty Clay Loam 

Pasture 

MFC 3 

51.0 

0 

36 

A2 

Silty Clay Loam 

Pasture 

Day 168 

97.17 

1.94 

36 

B1 

Silty Loam 

Soybeans J_ 

MFC 3 

52.0 

1.98 

180 

B1 

Silty Loam 

Soybeans 

Day 1 68 

87.14 

5.79 

154 

B2 

Silty Loam 

Smooth Bare 

MFC 3 

52.0 

1.89 

190 

B2 

Silty Loam 

Smooth Bare 

Day 168 

66.40 

5.50 

170 

Cl 

Si 1 ty Loam 

Soybeans _|J_ 

MFC 3 

57.0 

0.55 

114 

Cl 

Silty Loam 

Soybeans JJ_ 

Day 168 

66.33 

6.07 

102 

C2 

Sandy Loam 

Soybeans JJ_ 

MFC 3 

25.0 

0.0 

60 

C2 

Sandy Loam 

Soybeans 

Day 168 

44.90 

4.97 

60 


'it 

Field number corresponds to thos'’; noted on Fig. 2.9- 

'^MFC "N" — "N" refers to moisture condition in Table 1.3. 

Day "X" — "X" refers to Julian day in dynamic soil water-accounting model. 
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168. The means and standard deviations of the identified fieids are 


given in Tabie 2.1. Image graytone is proportional to moisture content 
in Figure 2.13; the upper image from the simple water budget is 
comparable (MpQ 10} to 30} drier) in mean moisture to the lower Image 
from Julian day 168 as produced by SWAM. The impact of the addition of 
within-fieid variance In moisture Is Immediately evident. The impact 
of crop type in SWAM Is exemplified by comparing Fieids B1 and B2 In 
both Images. In Figure 2.13a, both fields have a mean 0-3-cm soil 
moisture of 52} of field capacity while SWAM yields a mean moisture of 
87} and 66} for the soybean and bare fields respectively (with both 
fields having identical rain histories). The impact of soil type in 
both models Is most apparent when comparing portions of the Soybean 
Field C, which Is partly located on a silt loam (Cl) and a sandy loam 
(C2). 

The final moisture distributions for the radar image simulations 
detailed In Section 3 are shown In Figure 2.14. The dynamic soil water 
accounting model w?~ initialized on Julian day 153 and the output 
sampled every five days. The resultant distributions were then 
examined and the three most closely approximating moderately dry, 
moist, and wet soil surface conditions were selected for radar Image 
simulation. The distribution presented for day 173 In Figure 2.14 
shows that a minimum of approximately 4.0} of all grid cells contain 
soils at saturation and are potentially under standing-water 
conditions. This does not seem unreasonable when considering that a 
minimum of 5 cm of rain fel I across the data base on the previous day. 
In large part, the saturated conditions exist beneath the center of the 
storm track where rainfall approached 9 cm and for poorly drained clay 
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Figure 2.14. Actual soil moisture distributions as output by the 
dynamic soil-water accounting model on Julian days 
158, 168, and 173 for 20 m grid cells. 
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soils. The effect of this local flooding Is apparent on the simulated 
radar Imagery In Section 3.0. Image presentations of 0-3-cm soil 
moisture distribution for Julian days 136, 168, and 173 are shown on 
Figure 2,15. 
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a) day 1 58 


b) day 163 

Figure 2.15- Image display of 0-5 cm soil moisture distributions 

established by the dynamic soil water accounting model 
on Julian days: 158, l63, and 1 3. Percent of field 

capacity is proportional to image intensity. 
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c) day 173 
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3.0 RADAR SIMULATION STUDY 


The objective of this series of radar Image simulations Is to 
reexamine prior conclusions regarding soli moisture estimate accuracy 
Cl»23 with respect to the resolution of an. orbital C-band radar for the 
more complex spatial dlstrlbutlons-of soli moisture generated by SWAM. 
In addition the three SAR configurations previously modeled Cl»23 
with nominal resolutions of 20 m, 100 m, and 1 km (Cases 2, 4, and 8 
respectively. In Table 1.2), an orbital real -aperture radar (RAR) with 
a nominal resolution of 3 km (Case 12 in Table 1.2) is also simulated. 

The slmul t:tIon*s data base and procedures are Identical to those 
detailed In the previous study Cl»2J except for the terrain category 
reassignments and tha variability In moisture distribution Imposed by 
SWAM. The mean backscatter coefficient for each grid cell Is computed 
from empirically determined functions for each target class relating 
local Incidence angle, 0-5-cm percent of field capacity, surface 
roughness, and row direction to the mean radar backscatter coefficient. 
The effects of row direction on radar backscatter are defined as those 
that are experimental ly measured for dry- 1 and farming practices 
followed within the simulated region, and row directions are randomly 
distributed among fields with a 0.5 probability of being either 
parallel or perpendicular to the simulated radar look direction. 

The Julian dates of simulated satellite overpasses and the 
system-design options of the simulated Images are given In Table 3.1. 
The high-i-esol utlon SAR systems (20-m and 100-m resolutions) are 
sl.nulated for all three moisture dlstrlbuTl ans (Figures 2.14 and 2.15) 
In order to examine the effects of moisture level and relative moisture 
distribution function on soil molslure estimate accuracy. The effects 
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TABLE 3.t 

Simulated System Design Cases and Julian Dates of Satellite Overpasses 





of varying radar rasolutlon and the resolution of the resultant soli 
moisture estimate are examined for ail system design configurations 
using the moderately moist soil water distribution on Julian day 1S8 as 
a standard for comparison. 

The simulated Images for each condition In Table 3.1 are presented 
In Figures 3.1 to 3.9. The Image tone (brightness) Is linearly scaled 
to received power In dB; the log transform Is made In order to present 
the full dynamic rwtge of received pow«- such that point sources 
saturate the image and shadows are In the noise. 

3.1 jjfflcts of Changing Soil Moisture (>;ndltlQns 

This section will concentrate on an analysis of the Images 
presented In Figures 3.1 to 3.6. An examination of Figures 3.1 to 3.3 
from an orbital SAR with 20-ra resolution for Julian days 158, 168, and 
173 reveal several interesting features resulting from SWAM. 

The storm model produced convective rainfal I across the upper 
quarter of the Image swath on Julian day 154 and across the lower half 
of the Image swath on Julian day 156. Simulated radar observation of 
the data base on Julian day 158 (Figure 3.1) shows the path of the 
second storm cell as a generally brighter Image tone while the passage 
of the first storm cell Is no longer dlscernable on the Image after 
four days of evapotransp I ration. In the moist areas In the lower half 
of the Image, the field patterns re not r distinct as they are In the 
drier upper portion of the Image. This Is due In part to the fact that 
Interfield variance In soil moisture becomes greater with drying 
conditions since SWAM resnonds differentially to crop type and crop 
stage. Also, the effects of local slope on radar backscatter act to 
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Figure 3.1. Simulated radar image A1 for 20 m resolution on 
Julian day 158. 
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Figure 3.2. Simulated radar image A2 for 20 m resolution on 
Julian day 168. 
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Figure 3.3. Simulated radar image A3 for 20 m resolution 
Julian day 173. 
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Figure 3.4. Simulated radar image B1 for 100 m resolution 
on Julian day 158. 
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Figure 3.6. ^’nulated radar image B3 for 100 m resolution 
Julian day 173. 
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Figure 3.7. Simulated radar image Cl for 1 km resolution on 
Julian day 158. 
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reduce the discrimination of fle!d boundaries In the hilly lower half 
of the Image. This latter explanation Is also applicable to Figure 3.2 
for simulated Imagery for which 10 days of drying have caused modal 
soli moisture to be reduced by 20 percent of field capacity. 

On Julian day 172, a rather extensive storm has added 3 cm to 9 cm 
of rainfall over the test site, and simulated radar observation on the 
next day (Figure 3.3} shows that within the relatively flat region of 
the Kansas River floodplain (upper- left third of the Image) the radar 
response k; sensitive to conditions related to soil type and also is 
dependent upon field-control led conditions. The roughly triangular 
dark region In the left center of the Image is caused by the more 
specular reflection from poorly drained clay soils which are at 
saturation in the surface layer or are covered with standing water. 
This Is Interesting In that a "blind” moisture estimation algorithm 
(such as the one used Cl ,23) which requires only received power and 
range to estimate soli moisture will predict that this area Is dry soil 
when. In fact. It Is saturated or flooded soil. This effect has boen 
experimentally observed by truck-mounted scatterometers for these soils 
after rains of similar magnitude In 1977 C233* AppI Ication of the 
moisture prediction algorithm given by Eq. 3.1 to Figure 3.3 shows 
this to be the case. In Figure 3.9, 20-m resolution radar estimates of 
soil moisture distribution are plotted for Julian days 138, 168, and 
173. For Julian day 173, a comparison of actual and predicted moisture 
distributions. Figures 2.14 and 3.9 respectively, reveals that roughly 
15j of the 20-m x 20-m grid cells reflect serious underestln»tos of 
0-3-cm soil moisture. The clay loam soil with poor drainage occupies 
13$ of the 20-m x 20-m grid cells, and between Julian days 168 and 173 
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Figure 3.5. Distributions of 0-5 cn soil moisture estimated 
from 20 m resolution radar images (Case A) on 
Julian days 158, 168, and 173. 


69 







the number of grid cells classified as water doubles from 2.67X to 5,1% 
of the data base, reflecting the magnitude of local flooding within 
poorly drained fields. 

A 

The absolute difference between the estimated soli moisture MpQ 
predicted by Eq. 1.1 and the "actual" moisture MpQ as produced by SWAM 
on a given day Is defined as estimate error E. The cumulative 
distribution of E as calculated for each 20**m x 20-m grid cell element 
Is plotted versus the magnitude of E In Figures 3.10 to 3.13 for each 
simulated satellite overpass date. Thus, for a 20-m resolution SAR, 
Figure 3.10 shows that between 671^ and 1\% of ai I 20-m x 20-m grid cel 1 
elements (of 800,000 comparisons) have absolute estimate error E less 
than 30^ depending on the overpass date. In Figure 3.10, the lack of a 
significant difference between overpass dates (moisture conditions) for 
E < 40} can be attributed to the compound Impact of positioning errors 
(In matching ApQ to Mpr with a +/- 100-meter accuracy), the effects of 
local slope on radar moisture estimate, ond the mechanics of the blind 
comparison procedure wherein the error function Is derived from all 
800,000 comparisons regardless of land-use class. In sudi a 
comparison, all grid cells containing non-agr I cultural or non-forest 
categories will have a large E by definition, since actual moisture Is 
undefined for such categories (roads, buildings, water, etc.) and 
these kinds of categories comprise approximately 8} of the data base. 

Perhaps a more realistic evaluation of radar soil moisture 
estimate accuracy Is presented In Figure 3.11 for the same 20-m radar 
resolution and moisture conditions as those shown In Figure 3.10. In 
tills case, all non-agr leu I tural grid cell elements have been excluded 
from the calculation of tne error dHtrlb’-tlon function. In addition. 
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Figure 3.10. Soil moisture estimate accuracy for each overpass for 
radar with a 20 m resolution (Case A) when compared 
to mean 0-5 cm moisture present in 20 m x 20 m grid cells 
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Figure 3,11. Soil moisture estimate accuracy for agricultural 
categories on the floodplain for each overpass of 
a simulated radar with 20 m resolution (Case A) when 
compared to mean C-5 cm moisture present in 
20 m X 20 m grid cells. 
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the ONnparlson Is limited to the relatively flat floodplain where mean 
registration error between MpQ and Mpc Images Is observed to be +/- 20 
m. For the remaining 183»000 20-m x 20-m grid cell comparisons, 
between 75$ and 86$ have an absolute estimate error of E < 30$. When 
expressed In terms of gravimetric moisture, an estimate error E of 30$ 
corresponds to +/- 2.7$, +/- 7.3$, and +/- 9.2$ for loamy sand, silt 
loam, and silty clay loam, respectively (Table 3.2). 

Figures 3.12 and 3,13 show comparable results for the accuracy of 
soil moisture estimates using a 100-m resolution radar (System (^se 4). 
The general level of the cumulative function shows that roughly 10$ 
more 20-m x 20-m grid cells are correctly estimated for a given error 
threshold than can be correctly estimated from flie 20-m resolution 
radar (System Case 2). This result confirms that found for the 
simplified, low spatial frequency moisture conditions obtained before 
the addition of SWAM Cl»23 as shown In Figures 1.2 and 1,3. 

In addition, the highest estimate accuracies are obtained for the 
Intermediate moisture condition present on Julian day 158. This s 
most apparent in Figures 3.11 and 3.13 which consider only agricultural 
grid cells on the river floodplain. In botn figures, the difference In 
cumulative percent between Julian days 158 and 168 Is not statistically 
significant; however, estimate accuracy Is shown to be considerably 
less for day 177. This Is a consequence of the orevlously noted local 
saturation and par'Ial flooding of some regions o' the floodplain, with 
these areas apparently accounting for about 10$ of the total regl_>n. 
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TABLE 3.2 


Estimate Accuracy Levels: Comparison 
of Absolute Error in Percent of Field 
Capacity to Percent Gravimetric 
Moisture for Loamy Sand, Silt 
Loam; and Silty Clay Loam 


± Estimate Accuracy in Percent 


Gravimetric Moisture 

Fit-ld Capacity 

Loamy 

Silt 

Silty 


Sand 

Loam 

Clay Loam 

10% 

0.9 

2.4 

3.1 

20% 

1.8 

4.8 

6.2 

30% 

2.7 

7.3 

9.2 
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Figure 3.13. Soil moisture estimate accuracy for agricultural 

categories on the floodplain for each overpass of a 
simulated radar with 100 m resolution (Case B) when 
compared to mean 0-5 cm moisture present in 20 m x 
20 m grid cell s. 
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3.2 Imoa-zt of Radar Resolution and Estimate Resolution on Soil 

Moisture Estimate Accuracy 

The simulation procedure allows the prediction of an estimated 

A 

soil moisture MpQ for any desired grid-cell size from the radar Imagery 
(even If the estimate Is Invariant over distances less than the sensor 
resolution). The result Is two dimensional itios of resolution which 
can be evaluated with respect to estimate accuracy. The first 
considers only the effect of changing the radar resolution with respect 
to a fixed concept of the ”ground-truthed" data base (actual Mpc) while 
the second considers the effects of changing the resolution of the 
moisture estimate to conform to a redefined actual moisture 

distribution (as averaged Into larger effective grid cel is — 100 m x 
100 m. 1 km X 1 km, or 3 km x 3 km). 

A 

When the soil moisture estimated from simulated radar Imagery Mpg 
on Julian day 138 at radar resolutions of 20 m, 100 m, and 1 km Is 

compared to the actual moisture MpQ from SWAM within each of 800,000 

20-m X 20-ra grid cells, the resul isnt error distribution functions are 
found to be those shown In Figure 3.14. The estimate error Is shown to 
be the most tightly distributed around zero error for the 1-km 
resolution radar and most broadly distributed for the 20-m resolution 
radar. 

The shapes of the distributions In Figure 3.14 also exhibit 

Instructive trends. The finer resolution radars (20 m and 100 m) 
exhibit broad shoulders where soil moisture Is underestimated; 
examination of the position of these errors reveals that they are 
spatially related to the shadowed, backslopes of steep hillsides and to 
the presence of deciduous trees (no penetration of deciduous trees Is 
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assumed for the target-sensor Interaction model at 4.75 GHz In summer). 
These errors are not observed at the 1-km resolution, which reflects 
the fact that these target conditions (which result In a decrease of 
local ) are general ly much smaller than approximately 1 km^ In 
extent. On the other hand, the 1-km-resolutlon radar yields a small, 
but broad, secondary maximum In the error distribution for an 
overestimate of soli moisture. This peak Is related to the presence 
and distribution of point targets and cultural features, which behave 
In a fashion similar to corner reflectors and saturate the received 
signal. For a coarse resolution system, the Impact of these targets Is 
averaged over a correspondingly larger aroa and hence the estimate 
error attributed to such features Includes more pixel elements than are 
occupied by the actual feature* 

Thus, a reasonable conclusion is that adequate resolution for soil 
moisture sensing Is limited, at the fine resolution end, by the size 
and slope distribution functions of hills and the size and shape 
distributions of woodland; and limited at coarse resolutions by the 
distribution functions of hard point targets which have a large radar 
cross- sect i on. 

The absolute error (distance from zero error In Figure 3.14) In 
the sol l-molsture estimate uS computed from 20-m x 20-m grid cell 
comparisons Is shown for 20-m, 100-m, and 1-km radar resolutions in 
Figures 3.15 and 3.16 for the full data base containing 800,000 grid 
cells and the agricultural floodplain containing 183,000 grid ceils, 
rospectively. When compared to Figures 1.2 and 1.3, those figures show 
that the addition of within-fleld variance of soli moisture produces a 
negligible effect on the resultant error functions. The spatial 
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Absolute Difference Between Actual and Predicted Mpc 


Figure 3.15. Soil moisture estimate accuracy for a satellite 
overpass on Julian day 158 for radar resolutions 
of 20 m, 100 m, and 1 km when compared to mean 
0-5 cm moisture present in 20 m x 20 m grid cells. 
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Figure 3.16. Soil moisture estimate accuracy for a satellite 
overpass on Julian day 150 for radar resolutions 
of 2" m, 100 m, and 1 km wher compared to mean 
0-5 cm soil moisture present in 20 m x 20 m 
grid cells. 
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avm^aglng of the 100-m and 1-km resolution SAR Is sufficient to offset 
the Intact of the high** frequency spatial variability In MpQ (with a 
standard deviation of around a zero mean). Furthermore, the effect 
of within-fleld variance In soli moisture is not discernable at a radar 
resolution of 20 m, even when considering only the 183,000 grid-cell 
comparisons on the agricultural floodplain In Figures 1.3 and 3.16. 
This may, however, be a spurious result due to the relative positioning 
accuracy of Mpc to Mpc on the order of +/- 20 ra on the floodplain 
to V- 100 m In hilly areas. In addition, the effects of signal fading 
at a 20-m radar resolution with 12 looks lead to a greater estimate 
uncertainty than that due to the local variance In soil moisture. 

The above conclusions are also confirmed at a radar resolution of 
3 km for the RAR system ((Dase 12 In Table 1.2). The cumulative percent 
of the total data base as a function of absolute estimate error Is 
shown In Figure 3.17 from 100-m rr'^olutlon SAR, 1-km resolution SAR, 

A 

and 3-km resolution RAR; estimated moisture MpQ Is compared to MpQ on a 
100-m X 100-m grid cell basis. For the moderate soil moisture 
conditions prevailing on Julian day 133, the estimate accuracy of the 
3->'m system Is effectively equivalent to that of the 1-km resolution 
system. 

The effects of changing the grid-cell scale of the estimated soil 
moisture are shown In Figure 3.18. All results are for the 1-km 
resolution SAR on JulMn day 158, the only difference being the 
effective resolution at which soil moisture prediction Is being 
compared to the actual moisture distribution producer by SWAM, As the 
size of the comparison grid cell increases from 20 m x 20 m to 100 m x 
100 m, the net estimate accuracy Increases and the distribution becomes 
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Figure 3.17. Soil moisture estimate accuracy for a satellite overpass 
on Julian day 158 for radar resolutions of 100 m, 1 km, 
and 3 km when compared to mean 0-5 cm moisture present 
in iOO m x 100 m grid cells. 
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Estimate Error. Percent of Field Capacity 


Figure 3.18. Soil moisture estimate error resulting from 

comparing that moisture estima*''d by a 1 km resolution 
radar on Julian day 158 (Case Cl) to the mean 0-5 cm 
moisture present in 800,000 20 ir, x 20 m grid cells; 
in 32,000 10 m x 100 m grid cells; or in 325 1 km x 
1 km grid cells. 
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less skewed with respect to zero estimate error. Further change In the 
distribution as comparison grid-cell size Is Increased from 100 m x 100 
m to the sensor resolution size of 1 km x 1 km Is In the same direction 
but Is not significantly different In magnitude. The cumulative effect 
of changing the sol l-molsture estimate scale Is seen In Figure 3.19 for 
the same data as Is plotted in Figure 3.18. Thus for a 1-km resolution 
radar, the most accurate estimate Is obtained relative to the actual 
moisture distribution when estimated soil moisture Is considered as 
applying to the mean moisture existent within grid cells no smaller 
than 100 ra x 100 m. For the given spatial distributions of land cover 
classes In the data base, as representative of dr /-I and farming in 
east-central Kansas, the use of a 1-km resolution SAR to estimate 
moisture at a scale of 100 m x 100 m results In a loss In accuracy o-^ 
only 1$ to 35t, for a given acceptable error threshold as compared to 
using a 1-km x 1-km scale to estimate soil moisture. 

in a statistical sense, this result Is shown to be valid also for 
the 3-km RAR as shown in Figure 3.20. The reasons for this are not 
clear to the authors; however. It Is felt to be related to the scales 
of natural variability In the constituent scene components. 
Agricultural fields have dimensionalities on the order of hundreds of 
meters, the magnitude of local slope varies over lateral dimensions of 
tens to hundreds of meters In the data base area, wooded areas within 
the data base have lateral extents that range from tens to hundreds of 
meters, mappable soil-type typically varies at scales of hundreds to 
thousands of meters. in response, near-surface soil moisture Is 
observed to vary at all scales from millimeters to kilometers as 
controlled by the above conditions as well as microenvironmental 
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Percent of Data Base with Less Error 


ORIGINAL PAGE ^ 
OF POOR QUALITY. 



Absolute Difference Between Actual and Predicted Mpc 


Figure 3.19. Cumulative effect on moisture estimate accuracy of 
tt changing the size of prediction grid cell from 20 m x 

i| 20 m to 1 km X 1 km for radar resolution of 1 km 

I (Case C) on Julian day 158. 
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OfflGlNAl- PAGE B 
OF POOR QUALITY 



Figure 3.20. Effect on soil moisture estimate accuracy of changing 
the size of the prediction grid cell from 
100 m X 100 m to 3 km x 3 km for a radar resolution 
of 3 km on Julian day 150 (Case D). 
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condttfons such as surface roughness, shading, etc. (whose net effects 
are modeled as a random local variance with a stancard deviation of 6% 
of field capacity). 

3.3 Effects of Variance In Soil Moisture on Sei isQr_fiesolutJon 

For any given area or moisture condition, the optimum sensor 
resolution appears to be the one that most closely matches the ground 
scale at which sol I moisture becomes laterally decorreiated. When 
moisture is controlled largely by regional factors such as climatic 
events (rainfall) and soil-type distribution, the coarse resolution 
systems (1-km SAR and 3-km RAR) seem to provide optimal estimate 
accuracy. Within the closed system described by the simulation data 
base and according to the assumptions inherent in SWAM, this condition 
Is well approximated by the saturated-to- very-moist soil conditions 
that can be expected to prevail for several days after a significant 
rainfall event. As time progresses, more localized environmental 
factors such as surface slope, crop type, and stage of growth play an 
Increasingly significant role In controlling the spatial distribution 
of soil moisture; hence, the optimal sensor resolution for soil 
moisture mapping can be expected to be on the order of hundreds of 
meters . 

Indeed, this Is shown to be the case for tills series of radar 
image simulations. Figure 3.21 shows that for the moderately moist to 
wet moisture conditions (Mpc from 85% to 120$) present on Julian day 
158, the actual moisture contritions are most accurately reproduced by 
the coarsest resolution systems (1-km SAR and 3-km RAR) at the high end 
of the moisture range and by the 100-m resolution SAR over the median 
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Figure 3.21. Comparison of actual percent of field capacity on Julian day 158 with that 
estimated from radar with resolutions of 20 m (Case Al), 100 m (Case Bl), 



moisture range. The 20-tieresol utlon SAR Is far too sensitive to local 
target-sensor Interactions which Introduce "noise” components Into the 
resultant soil moisture estlmate^and when combined with procedural 
problems In accurately registering the moisture estimate relative to 
the Input MpQ distribution, yields thfe relatively lower estimate 
accuracy observed In Figures 3.14 to 3.16. 

However, regardless of radar resolution, the impact of local 
withln**f leld variability In soil moisture seems to have a negligible 
Impact on soil moisture estimate accuracy. This apparently results 
from the fact that 

•local within-fleid variance In soil moisture Is normally 
distributed with a zero mean, 

•Rayleigh fading Is a larger source of local (20-m x 20-m} 
soil moisture estimate error for a small number of looks 
(N < 23) , and 

•other local effects (such as the presence of hills, 
canopy cover types, and row direction effects) 
are more significant sources of soil moisture 
estimate error from radar Imagery. 

4.0 CPNCLUSIPMS 

The major conclusions of this study may be summarized as: 

1) Local near^surface sol l-molsture content was estimated with the 
greatest accuracy from radar Imagery produced by systems with 
resolutions coarser than 100 m. Of the three such system 
configurations tested (100-m resolution SAR, 1-km resolution SAR, and 
3-km resolution RAR) the optimal resolution was found to be linked to 
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the Interaction between general soli moisture level t^nd the spatial 
distributions of environmental factors controlling soil moisture 
dynamics. Thus« the coarsest resolutions (1 km and 3 km) performed 
best for very moist soil conditions where surface moisture content is 
primarily deten !' jd by soil type and drainage characteristics, while a 
sensor resolution of 100 m displayed equivalent or better moisture 
estimate accuracy during drying conditions as differential dry-down 
rates of various crop canopies Introduced significant interfleld 
variability Into the soll-molsture dIstrloutJon. 

2) The effects on soil moisture estimate accuracy of 
high-frequency spatial variance In true soli moisture are negligible 
for radar resolutions of 100 m or larger If the variance is randomly 
distributed with a zero mean. For a radar resolution of 20 m with 12 
looks, the uncertainty In moisture estimate as related to fading 
confidence Interval Is greater than that due to local target 
variability modeled as a Gaussian with a standard deviation of +/-6 
percent of field capacity. 

3) The upper limit of desirable resolution should be bounded by 
considerations of field-size distribution for a given agronomic region, 
the soil moisture range of Interest for a specific application, and the 
spatial density of non-dlstrlbuted targets with a large radar 
cross-section. Thus, the radar resolution cell should be less than the 
modal agricultural field size of interest, and small enough to allow 
the discrimination of point targets as sources of soll-molsture 
estimate error (especially for dry soll-molsture conditions). 

The significance of the above conclusions Is, of course, limited 
by the authenticity of the closed system described by the dynamic 
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envlrorenontal model, the static terrain model, the target-sensor 
Interaction model, the radar Image formation process, and the 
evaluation procedures applied to the resultant Images. The static 
terrain model Is certainly appropriate for the simulated area, except 
that the Inclusion of certain types of scattered point targets Is 
somewhat arbitrary at a 20-m x 20-m grid-cell size. The dynamic soil 
water accounting model, while not rigorous, correctly predicts the 
direction and gross magnitude of changes In near-surface soil moisture. 
Tho radar backscatter models have been found to be statistically 
accurate descriptors of mean conditions; the explicit addition of crop 
growth stage to backscatter I ng models Is a planned Improvement that 
will require more detailed empirical evidence as well as a better 
theoretical understanding of the Impact of crop phenology on scattering 
and absoprtlon. The Image formation models have been found to be quite 
good for side- looking airborne and orbital radar, which Incorporate 
Rayleigh fading statistics for Images containing a small number of 
Independent samples C26,27]]. 

Thus, the results are expected to represent fairly the conditions 
applied. It Is, however, recommended that the geographic extensibility 
of these conclusions be tested for a larger or geographically 
dissimilar data-base, especially with respect to these conditions. 

•agricultural field-size distributions with mean field- 
size larger and/or smaller than that modeled, 

•more Intricate spatial distribution of agricultural fields 
and woodlands such as those characteristic of agronomic 
areas In the more humid eastern half of the United States 
and Europe, 

•denser spatial distributions of cultural features and point 
targets characteristic of more densely populated yet 
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agriculturally significant regions. 

Also, given the apparent adequacy of coarse- resolution radar 
systems with low power requirements, consideration should be given to 
modeling both a cross-polarized system and a combined rodrr/radlome-*’er 
system for purposes of estimating soil moisture. 
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